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Abstract 

We evaluate the average and the standard deviation of the supercurrent in superconducting 
nanobridges, as functions of the temperature and the phase difference, in an equihbrium situation. 
We also evaluate the autocorrelation of the supercurrent as a function of the elapsed time. The 
behavior of supercurrent fluctuations is qualitatively different from from that of the normal current: 
they depend on the phase difference, have a different temperature dependence, and for appropriate 
range their standard deviation is independent of the probing time. We considered two radically 
different filaments and obtained very similar results for both. Fluctuations of the supercurrent can 
in principle be measured. 



PACS numbers: 74.40.-n, 74.78.Na, 72.70.+m 



I. INTRODUCTION 



Josephson junctions are of common use in many technological devices and their behavior 
as a circuit element is described in textbooks;-"- the influence of thermal fluctuations on the 
normal current is described as Johnson noise. The voltage due to thermal noise was studied 
by Ambegaokar and Halperin.- 

In this study we are interested in the equilibrium fluctuations of the supercurrent. Since 
supercurrent is actually an equilibrium variable rather than a diffusion process, we may 
expect — and will indeed flnd — that its fluctuations are qualitatively different from those of 
the normal current. Early works on fluctuations in junctions^ state that there is no noise in 
supercurrent to quadratic order in the temperature. On the other hand, Averin and Iman>^ 
found that supercurrent fluctuations in mesoscopic contacts are large on the scale of the 
classical shot noise. 

The kind of junction considered in this article is a fllament close to the critical tempera- 
ture, which can be described by means of the Ginzburg-Landau model. If the entire fllament 
has a critical temperature that is above the experimental temperature T, the fllament may 
be regarded as a constriction; if parts of the fllament have a critical temperature above 
T, it may be regarded as an SNS junction. In the case that fluctuations are ignored, the 
junction-like behavior of a fllament has been studied for static'' and for dynamic^ situations. 

The situation we will study here is a case of dynamic equilibrium. We will consider a 
superconducting fllament that bridges between two "banks." At the banks fluctuations are 
negligible and the order parameter will have flxed equilibrium values, whereas along the 
fllament the order parameter and the electromagnetic potential fluctuate. In the absence of 
fluctuations, the current along the fllament would be given by the current-phase relation. 

II. METHOD AND DEFINITIONS 

We will use the time- dependent Ginzburg-Landau (TDGL) equations with Langevin 
terms, which will be handled numerically by means of flnite differences. We have described 
this method in detail in the past^ii^ and shown good agreement with statistical mechanics 
and with experiments. 

For a ID fllament we deflne the gauge-invariant order parameter ipi^) = 




2 



exp[(27ri/$o) Jq A{s')ds']4'{s), where s is the arc length, the "canonical" order parameter, 
A the tangential component of the vector electromagnetic potential and $o the quantum of 
flux. The ID-TDGL equation can be written as^ 



dip 



~ 2 h"^ d'^ 1 dw d 
2m ds"^ w ds ds 



^, (1) 



where m is the mass of a Cooper pair, w{s) the cross section of the filament, and u, a and 
P are material parameters; the sign of a determines whether the local critical temperature 
is above or below T. 

We will take the boundary conditions ^/'(0,t) = a/— a(0)//3(0), 'ijj{L,t) = 



a(X)75(Lyexp(z7), where L is the length of the filament; 7 is the gauge- invariant phase 
difference. 

The scaled values of a and (3 will be kept fixed in the present study. Since a and (3 are 
functions of temperature, the physical meaning of this exposition is that for each temperature 
the experiment is performed on a different sample. 

The normal current I]sr{s,t) and the supercurrent Is{s,t) are not separately con- 
stant along the filament. We define the supercurrent as the weighted average /c (t) = 



10 



m 



Is{s,t)w{s)~^ds/ w{s)~^ds. w{s)~^ appears as a natural weight in Refs. [9 
this work we consider uniform cross section only, so that this weight and the last term in 
Eq. ([1]) can be ignored. 



III. RESULTS 

We have examined two cases. One of them is a uniform filament with a = —h'^/mL'^] in 
the second case we took a{x) = —{h'^/'mL'^) cos{2itx/L), so that the middle of this filament 
is nominally normal. In both cases we took /3 = h'^w/rriL and the resistance of the filament 
as h/Aue^. The filament was divided into 30 computational cells and evolution was followed 
in steps of duration 1.3 x 10~'^umL'^ / h. The first 10^ steps had the purpose of relaxation 
to typical values of ip and then 12 x 10^ steps were used for averaging. Our results do 
not depend appreciably on the resistance or the steps duration. The averages shown in the 
figures are averages over time and Iso denotes the supercurrent in the absence of thermal 
fluctuations. 

Using BCS, dirty limit^ and free electron gas approximations, our choices for a and /3 
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FIG. 1: Average deviation of the supercurrent from the value that would be obtained without 
thermal fluctuations, as a function of the phase difference and for several temperatures (solid 
curves). The temperature is O.ln'^h? /mL^kB, where n is the number marked next to each curve. 
For comparison, we have also drawn the curve — /50/3 (dashed line). 

can be inverted and regarded as choices for the geometric parameters. We obtain 

L2 = Tih'^kFie/GrnkBiTc - T) , w = 0.51L/nJl , (2) 

where kp is the Fermi wavevector, ig the mean free path, the critical temperature of the 
banks and ng the electron density. 

A. Uniform filament 

Figure [T] shows the average deviation of Is from J50 as a function of 7 for several tem- 
peratures. We see that fluctuations not only lead to variance of the supercurrent, as in the 
case of normal current, but also lead to a shift of the average value. Is and I so are negative 
for < 7 < vr, so that \ {Is) \ < l-^sol- For a material with critical temperature of the order of 
IK, the highest temperature in Fig. [1] corresponds to a filament length of the order of 30 nm 
and the highest current deviation, to the order of 10 nA. The 7-dependence of (Is) —Iso has 
some resemblance with that of Iso, which has been included in the figure for comparison. 
Figure E] presents (Is) — Iso as a function of temperature; it is apparent that this deviation 
behaves differently for different values of 7. 

Figure [3] shows the standard deviation of Is as a function of 7. Note that whereas 
the standard deviation of Johnson noise is inversely proportional to the square root of the 
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FIG. 2: Deviation of the supercurrent from the fluctuation-free value, as a function of the temper- 
ature and for several phase differences, marked next to each curve. 
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FIG. 3: Standard deviation of the supercurrent, as a function of the phase difference, for several 
temperatures. The temperatures are the same as in Fig. [TJ 

probing time (provided that it is long compared to h/ksT), the standard deviation of Is 
is independent of the probing time (provided that it is short compared to the decoherence 
time that we will find below). In addition, the standard deviation of Is depends on the 
phase difference, whereas that of In does not. On the other hand, the present result does 
not support the scenario assumed in Ref. il2i, according to which fluctuations of the order 
parameter just scale the supercurrent, while the shape of the current-phase relation remains 
fixed; if this were the case, the standard deviation of Is would be proportional to Iso- 
Figure H] shows the temperature dependence of the standard deviation of Is- Whereas 
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FIG. 4: Temperature dependence of the standard deviation of the supercurrent. The lowest curve 
is for 7 = 0, the highest for 7 = vr, and the curves in between are in steps of 7r/6. 
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FIG. 5: Autocorrelation function of 1$, as a function of the elapsed time t', for temperature 
2.blT? /niLp'kB and 7 = nvr/S, n = 0, 1, 2, 3. The other parameters are as in Fig. [TJ Inset: spectral 
density for 7 = 0, 7r/2, tt. 

for the normal current the standard deviation is proportional to T^/^, for the supercurrent 
this scaling occurs only in the case 7 = tt; as 7 decreases towards 0, the initial slope of the 
curves decreases. 

Figure O shows the autocorrelation of the supercurrent, K(t') = {{Is(t + t')Is(t)) " 
(/5)^)/((/|) — (Is)"^) for < 7 < TT. We see that, the smaller the value of 7, the shorter 
the typical time required to "forget" a previous value of Is- The inset shows the spectral 
density J{uj) = (I/tt) K{t) cos{ijjt)dt. 
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FIG. 6: Influence of temperature on the autocorrelation function. The dashed lines are for tem- 
perature close to zero and the continuous lines for 2.5h? /mL'^kB- The lower (blue online) lines are 
for 7 = 0, the middle lines for 7 = tt/2, and the upper (red online) lines for 7 = tt. 

If Is is measured many times during probing periods of length r and K{t) ~ 1, then 
each measurement can essentially be regarded as instantaneous and the standard deviation 
of Is is given by Fig. [3l with practically no r-dependence; if K{t) ~ 0, then fluctuations of 
Is essentially become white noise and the standard deviation of Is should decrease as r"^/^. 
According to Fig. |5l the crossover value of r (the "decoherence time" ) increases with 7 and 
is of the order of O.lumL'^/h (for L ~ lO^^m, this is of the order of lO^^'^s). 

The faster relaxation of fluctuations for small 7 is counterintuitive, since according to 
Fig. |3] the influence of thermal agitation is stronger for large 7, and we might expect this 
agitation to destroy any particular conflguration of the order parameter that results in 
a particular supercurrent at a given time. Figure [6] shows that the 7-dependence of the 
relaxation times of Is is not dominated by thermal agitation, but rather by the fluct nationless 
dynamics. At T ~ memory loss of particular conflgurations is faster for small 7; thermal 
agitation moderates the 7-dependence. It is interesting to note that in the case 7 = 
thermal agitation leads to delay of the relaxation. 

The results in Figs. [3] and |6] can be interpreted as follows. In the case 7 = the order 
parameter is pinned with the same phase at both boundaries, giving a large energy advantage 
to a uniform order parameter all along the filament. Therefore, the system is comparatively 
rigid: deviations from the equilibrium configuration are small and return to equilibrium is 
fast. In the case 7 = tt the order parameter is pinned with opposite phases at the boundaries. 
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FIG. 7: Reduction of average supercurrent due to thermal fluctuations in the case a = 
— (/i^/mL^) cos(27rx/L). Annotations as in Fig. [TJ For comparison, the curves for a = —t^jral?, 
T = OAh'^/mL^kB and T = im^/mL'^kB have been included as dotted lines. 

so that the system is frustrated and comparatively indifferent; as a consequence, deviations 
from equihbrium are large and return to it is slow. 

B. SNS junction 

Here we report the results for the case a = — (/i^/mL^) cos(27rx/L). Although taking 
a value of a that vanishes on the average may look as a drastic change, the results are 
remarkably similar to those of the previous section. 

Figure [7] shows the average deviation of Isi'j) from /so (7) for the same temperatures 
as in Fig. [TJ We note that in the present case the reduction of Is is smaller than in the 
case of a uniform filament. This behavior makes sense, since when the entire filament is 
superconducting fluctuations may only be expected to destroy superconductivity; on the 
other hand, when half of the filament is normal and supercurrent can only be expected 
due to proximity or fluctuations, fluctuations also have a supportive effect. For 7 = 7? the 
average current has to vanish by symmetry; the scattering in this value serves as a measure 
of the accuracy of our results. 

Figure [8] shows the standard deviation of Is{l) and compares it with that of a uniform 
filament. For small values of 7 the standard deviation is similar to that of Fig. [3l but 
as 7 approaches tt the standard deviation in the present case is noticeable smaller than 
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FIG. 8: Standard deviation of the supercurrent in the case a = —{h^/mL?') cos(27rx/L), for several 
temperatures. The temperatures are the same as in Fig. [TJ For comparison, the curves of Fig. [3] 
are redrawn as dotted hues. 
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FIG. 9: Temperature dependence of the standard deviation of the supercurrent. Sohd Hues: a = 
— {h? /mL'^) cos{2iTX / L); dotted: a = —h?/mL'^, redrawn from Fig. HI 

that of a uniform filament. This result may seem surprising, since in the present case 
superconductivity is more fragile and we might expect to larger fluctuations. A possible 
explanation could be that when a and /3 are both positive, large fluctuations of the order 
parameter in the middle of the filament are inhibited, leading to smaller fluctuations of Is- 
Figure [9] presents these results as functions of the temperature. In the present case the 
standard deviation does not become proportional to T^/^ even for 7 = vr. 

Figure [TH] compares the autocorrelation functions for the cases a{x) = 
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FIG. 10: Autocorrelation of the supercurrent, as a function of the elapsed time, for a = 

— [h"^ /mL'^) cos {2ttx/ L) and several values of 7. The other parameters are as in Fig. [5l The 
dotted lines are for a = —h'^/niL'^. 

— (/i^/mL^) cos(27rx/L) and a = —fv^/mL?. Again, we find that the difference is re- 
markably small, especially for small 7. Weakening of superconductivity leads to a slightly 
faster decoherence. 

IV. DISCUSSION 

We have evaluated numerically the thermal fluctuations of the supercurrent along fila- 
ments that bridge between two banks (superconducting pieces with dimensions such that 
fluctuations in them are negligible) in an equilibrium situation. One case we considered was 
that of a uniform filament, made of the same material as the banks; the second case was 
that of a filament equal to the first at the contact points, but normal in the middle. 

These fluctuations have non trivial properties, qualitatively different from those of normal 
current. On the other hand, in spite of the blatant difference between both considered 
filaments, the difference between the fluctuations in them is minor, suggesting that the 
results we have found are generic. 

Experimentally, a phase difference may be applied by connecting the banks so that to- 
gether with the fllament they become a closed circuit that encloses a known magnetic flux; 
the current can then be sensed through the fleld that it induces. Measurement of J50 is 
complicated by the fact that changing the temperature would also change the value of a. 
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Let us now discuss measurement of the standard deviation of Is- The probing time should 
be of the order of mL'^/h, which in view of Eq. (|2]) is of the order of kkpie/kBiTc — T), and 
ought to be accessible for temperature sufficiently close to T^. For this probing time, using 
our results in Fig. HJorlHl Eq. and the free electron gas expression for the resistivity, we 
estimate that the ratio between the standard deviation of Is and that of Johnson noise is of 
the order of (iekpY^'^'-i since this is not a small number, fluctuations of the supercurrent are 
not obscured by those of the normal current and could in principle be observed. 

The present study may be regarded as a feasibility test for the influence of supercurrent 
fluctuations. Analytic treatments that uncover the scalings and asymptotic relations, as well 
as the non equilibrium behavior, are still required. 
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Appendix A: Evaluation of {Is) using a transfer operator technique 



The equilibrium average of the supercurrent is given by 



(Al) 



$0 d-i 

where Z is the partition function, which has to be derived from the Ginzburg-Landau free 
energy 



wds[a\i)\^ + {l3/2)\i)\^ + {h^/2m)\di)/ds\ 



(A2) 



For a uniform filament and a < 0, writing s = Lt and (Rei/j, Imijj) = y/— a;//3 r, F becomes 



dt 



L(5m 



dr 

'dt 



+ V 



with 



V = {aL'm/h'){r' -ry2) 



(A3) 



(A4) 



Following Ref. Il3|, a function r{t) is interpreted as a microstate of the system and F as the 
energy of the system for that microstate. It follows that the partition function is 



Z = C Vrexp{-F/kBT) 



(A5) 
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where J Vr denotes integration over all functions r(t) with boundary conditions r(0) = ao = 
(1, 0) and r(l) = ai = (cos 7, sin 7), and C is an irrelevant multiplicative constant. Dividing 
the integral in Eq. (IA3p into N segments, N ^ 1, and introducing into Eq. (lASp we can 
write in this limit 

N-l 



Z = (^^^ j dri... dr^Siri - ao)(5(rjv - ai) exp[- '^^^^^^ . . . exp[ 



/(r2,ri)^ 



NS ' ' 
(A6) 

with /(r.j+i,r.j) = ^(rj+i) + (A^^/2)|rj+i — r^p, S = —kBTLf3m/wah'^, and the prefactor has 
been chosen for convenience. 

Noting that 6{r — r') = J2n ^n('^)^n(j'')! where {"^n} is any complete set of normalized 
eigenstates, Z becomes 

Z = ^^„,(ai)*:(ao) J rfr;v^:,M--- J dr,,(Ar/27r5) exp[-/(r,+i, r,,)/iV5] ■ • • 

n,n' 

X J dri{N/27rS)exp[-f{r2,ri)/NS]^4ri) . (A7) 
If the are chosen so that they obey the eigenvalue equation 

J dv,{N/2TiS) exp[-/(r,+i, r,)/Ar^]^„(r,) = exp[-e„/Ar^]^„(r,+i) , (A8) 

then 

Z = Y, ^n(ai)^:(ao) exp(-e„/5) . (A9) 

n 

Expanding \l'„(rj) in powers of rj — Tj+i around \I^„(rj+i), which is equivalent to an 
expansion in powers of N'^^"^, the integral in Eq. (lASP can be performed and is found 
to obey the eigenvalue equation^^ 

[-{Sy2)V^ + V]^n = en^n, (AlO) 

where is the Laplacian with respect to r. 

In polar coordinates r = r(cos6', sin 6'), the angular momentum operator is = —id/dO. 
Noting thair^ commutes with the "Hamiltonian" [— (5*^/2)^7^ + y], the can be chosen 
so that they are also eigenstates of L^, i.e. they can have the form \l/„ £[r(cos^, sin^)] = 
-Rn/('") exp(i£6'), where Rn/ is a real function and £ an integer. Introducing this form into 
Eq. (1A9P we obtain 

Z = 5^exp(2£7)Z, , (All) 
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with 

z^ = Y,RlA^)^M-^n,i/s), (A12) 

n 

where summation in Eq. (lAlip is made over all integers and in Eq. (lA12p over all the states 
with total angular momentum ^. 

Since the Hamiltonian is symmetric under the transformation r — )■ — r, we can also write 

oo 

Z = Zo + 2 ^ cos(£7)Z^ . (A13) 



1=1 



Applying Eq. (lAip we obtain 

lTa\ = -- 



(/^) = -^^^sin(£7KZ,. (A14) 



1=1 

In order to complete the evaluation, Eqs. (1A14P and (]A12p have to be supple- 
mented with the values of and Rn/iX). Since only for small values of e„/ there 
is an appreciable contribution, this is accomplished as follows. We start from a basis 
Hamiltonian Hb = — (5*^/2)^7^ + (A;^/2)r^, which has known eigenf unctions R^^{r) = 
Cnir\^\e-^''^l'^^ iFi{-n, |£| + l,kr'^/S) and eigenvalues ef^^ = Sk{2n + |£| + 1), where Cn/ 
is the normalization constant and i-Fi is Kummer's hypergeometric function, k is still a 
free parameter. We select a subspace of the Hilbert space which has a moderate number 
of low-energy eigenstates of Hb as a basis, project the true Hamiltonian into this subspace, 
and then diagonalize this truncated Hamiltonian. 

The value of k is chosen as follows. If the Hamiltonian were not truncated, its lowest 
eigenvalue would be independent of /c; for the truncated Hamiltonian and for a given value £, 
the lowest eigenvalue is effectively independent of k within a certain range and rises beyond 
this range. We choose k{€) in the middle of this range. 

Figure [H] compares between the computational method taken from Refs. and the 
method developed in this appendix. The comparison requires extrapolation to the continuum 
limit, — )■ oo. The convergence of {Is) to the continuum limit is rather slow, but the 
variance of 1$ seems to saturate for A^ ~ 20. For temperatures lower than OAh"^ /mL'^kB 
convergence in Eq. (1A14|) becomes slow. 
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FIG. 11: Supercurrent as a function of the length of the computational segments. The dashed lines 
extrapolate to the values obtained by the method described in this appendix. The temperature is 
2.bfi? /mLp'kB, the phase difference is marked next to each line and the other parameters are as in 
Fig. [TJ We kept 12 terms in expansion ()A14p . 
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